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LONG-TERM  GOALS 

The  long-term  goal  of  the  research  is  to  increase  the  physical  understanding  of  acoustic  propagation  in 
continental  shelf  and  slope  environments  in  the  50-4000  Hz  band.  This  includes  both  the  physics  of 
the  seabed  and  the  coupling  to  physical  mechanisms  in  the  water  column  in  complex  range-  and 
azimuth-dependent  littoral  waveguides. 

OBJECTIVES 

There  were  two  main  objectives  of  the  current  research.  The  first  objective  was  to  develop  and  test  a 
new  coupled-mode  approach  that  can  serve  as  the  basis  for  an  advanced  seabed  attenuation  analysis  of 
broadband  propagation  data.  The  second  objective  was  to  apply  the  coupled  mode  approach  to 
investigate  the  physics  of  scattering  from  a  rough  seabed  surface. 

APPROACH 

A  MATLAB  numerical  algorithm1  was  constructed  for  a  two-way  integral  equation  coupled  mode 
approach  that  was  originally  programmed  in  FORTRAN.2'3  The  MATLAB  algorithm  includes  a  more 
efficient  numerical  method  to  compute  the  range-dependent  mode  coupling  matrices.4'5  Then,  the 
MATLAB  algorithm  was  modified  to  explicitly  compute  the  forward  and  backward  going  coupled 
amplitudes,  as  discussed  in  Ref.  6.  Next,  the  modal  power  flow  within  the  trapped  spectrum,  into  the 
continuum,  and  into  the  backward  propagating  modal  spectrum  was  examined  using  the  master 
equation  from  statistical  mechanics,  as  discussed  by  Dozier  and  Tappert7  and  Colosi  and  Morozov.8 

The  propagation  and  scattering  theory  parallels  Ref.  6.  The  2-D  coupled  mode  equations  of  dimension 
2 N  for  the  forward  and  backward  modal  amplitudes  ( R  (r)  =  \  R’v  R™,—,  R°°N_V  R'N  \ 

and  R°(r) = R0^,..., ^+JV_P  in  cylindrical  coordinates  are 
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(i) 


pn  is  the  modal  source  function  for  the  nth  mode,  and  C  i s  the  mode  coupling  matrix  operator.  N  is 
the  total  number  of  modes,  including  both  trapped  and  leaky  modes.  G°  and  G"  are  the  Green’s 
functions  that  satisfy 


or  r  or 


S(r—r') 


where  kn(r )  is  the  horizontal  wavenumber  eigenvalue  for  the  nth  mode  at  range  r. 
neglected  in  Eq.  1,  the  result  is  a  one-way  equation  for  R°°. 


(2) 

When  R°  is 


Having  a  solution  to  Eq.  1  permits  the  modal  power  flow  to  be  calculated.  Following  Dozier  and 
Tappert7,  one  can,  using  the  Bom  approximation  and  neglecting  backscatter,  write  a  coupled  diffusion 
equation  (called  the  master  equation)  for  the  modal  amplitudes 


d/dr<Rn(0\2>=  ZC™II R10  f  ~ I  «,(r)  |2]- 2a  |  RJr) \\  (3) 

m 

where  an  attenuation  term  has  been  added,  as  discussed  by  Colosi  and  Morozov.8  We  will  show  the 
results  of  constructing  d  I  dr  <)  Rn(r)\2>  (LHS  of  Eq.  3)  and  X  CJ I  Rm(r)  \2 -\  R^r)]1  \-2a\  R^r)]1 

m 

(RHS  of  Eq.  3)  from  both  the  full  two-way  solution  and  the  approximate  one-way  solution  to  Eq.  1 .  If 
the  RHS  is  a  good  approximation  to  the  LHS  of  Eq.  3,  then  one  might  expect  to  observe,  in  the  limit  of 
small  attenuation,  an  equal  partitioning  of  modes  from  the  scattering.7 

WORK  COMPLETED 

The  work  completed  in  FY12  included  the  application  of  the  split  integral  equation  coupled  mode 
methodology  to  quantify  the  modal  power  flow  for  selected  seabed  roughness  realizations.  A  potential 
important  discovery  was  that  for  a  seabed  roughness  characterized  by  a  relatively  small  correlation 
length  the  energy  lost  from  the  trapped  modes  to  the  continuum  and  backscattering  can  be  less  than  the 
energy  saved  by  a  transfer  of  energy  within  the  trapped  spectrum  to  the  lower  order  modes  with  smaller 
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attenuation.  The  nature  of  this  type  of  seabed  scattering  can  cause  the  modal  energy  spectrum  not  to 
be  evenly  distributed  as  would  be  predicted  by  a  master  equation  approach  from  statistical  physics. 


RESULTS 


Selected  numerical  computations  from  the  IECM  approach  are  presented  to  illustrate  scattering 
mechanisms  resulting  from  a  rough  seabed.  The  background  waveguide  is  50  m  deep  with  a  stratified 
water  sound  speed  profile  (SSP)  over  an  infinite  half-space.  The  SSP  has  a  thermocline  depth  at  about 

'J 

20  m.  The  geoacoustic  parameter  values  describing  the  bottom  half-space  are  1645  m/s,  2.0  g/cm  ,  and 
0.4  dB/A.  for  the  sediment  compressional  sound  speed,  the  density,  and  the  attenuation,  respectively. 
The  maximum  (computational)  range  in  the  waveguide  is  10  km,  and  the  source  is  located  at  r  =  0  m 
and  depth  10  m.  The  roughness  wavenumber  power  spectrum  employed  to  generate  a  1-D  rough 
seabed  realization  is 


P{k)  = 


KLh2 

n{K2L+k2yn  ’ 


(4) 


where  h,  KL ,  y,  and  k  are  the  rms  roughness  height,  the  roll-off  wavenumber  given  by  the  reciprocal  of 
the  roughness  correlation  length,  the  scattering  exponent,  and  the  spatial  wavenumber,  respectively. 
For  the  current  study  scattering  effects  were  computed  at  y=2.0  and  h  =  0.5  m  with  Ki  varying  from 
0.0025  1/m  to  0.1  1/m.  For  the  numerical  computations  the  roughness  power  spectrum  is  employed  to 
generate  a  rough  seabed  for  3.0  km  <  r  <  7.0  km,  and  in  the  regions  0.0  km  <  r  <  3.0  km  and  7.0  km  < 
r  <  10.0  km  the  waveguide  is  range-invariant. 

Figure  1  shows  |  Rn  (r)  |  (adjusted  for  cylindrical  spreading)  obtained  from  the  full  solution  of  the  two- 

way  coupled  integral  equations  for  an  acoustic  frequency  of  200  Hz.  For  this  waveguide  and  frequency 
there  are  5  trapped  modes  and  8  continuum  (leaky)  modes  for  a  total  of  13  modes.  The  energy  in  the 
trapped  modes  is  clearly  shown  to  scatter  into  both  the  continuum  modes  for  3.0  km  <  r  <  7.0  km  and 
into  the  backward  propagating  modes  (modes  14-26  on  the  vertical  axis). 

Figure  2  shows  A (K£)  where  A  is  defined  as  the  2-way  depth-integrated  transmission  loss  (DITF)  minus 
the  adiabatic  (no  mode  coupling)  DITF  at  10  km.  A  becomes  negative  at  about  KL  =  0.06  1/m  and 
continues  to  decrease  for  increasing  KL.  An  examination  of  the  modal  amplitudes  showed  that  more 
energy  was  being  scattered  into  mode  3  than  into  the  continuum  and  into  backscatter.  The  distance 
between  interference  peaks  for  modes  3  and  4  is  ID34  =  2jt/|Re(k3)  -  Re(k4)|  «70  m.  For  small  KL,  ID34 
is  small  compared  to  l/KL.  As  KL  increases  the  spatial  correlation  length  decreases  and  becomes  on  the 
order  of  ID34  and  the  interference  length  for  other  modal  combinations  within  the  trapped  spectrum. 
More  generally,  when  the  modal  interference  distance  is  small  compared  to  l/KL  the  large  number  of 
modal  interference  cycles  over  a  correlation  length  allow  for  a  significant  transfer  of  energy  between 
these  modes  to  occur.  Such  scattering  will  be  favored  over  pairs  of  modes  that  have  a  larger  ID  value. 
However,  as  KL  increases  the  correlation  distance  decreases  and  the  preference  for  scattering  into  the 
continuum  from  the  trapped  spectrum  diminishes.  To  establish  the  uncertainty  in  Fig.  2,  the  coupled 
mode  equations  were  solved  for  20  roughness  realizations  for  selected  KL  values.  The  result  was  that  a 
=  0.76  at^i  =  0.1  1/m  and  a  =  0.48  at  KL  =0.018  1/m. 
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Figure  3a  shows  dldr  <|  7^(r)|2>,  where  the  modal  amplitudes  were  those  used  to  construct  Fig.  1. 

Figure  3b  shows  dldr  <|  7^(r)|2>  computed  with  the  RHS  of  Eq.  3.  A  comparison  of  Figs.  3a  and  3b 

suggests  that  the  RHS  of  Eq.  3  overestimates  the  amount  of  acoustic  intensity  radiated  into  the 
continuum;  however,  the  RHS  of  Eq.  3  underestimates  the  amount  of  energy  transferred  into 
backscatter.  Some  of  the  differences  between  Figs.  3a  and  3b  may  be  ascribed  to  the  inclusion  of 
backward  propagation  into  the  RHS  of  Eq.  3,  which  is  neglected  in  the  original  derivation  of  Eq.  3  by 
Dozier  and  Tappert.  Figure  3b  suggests  an  equal  partitioning  of  modes  (trapped  and  continuum) 
whereas  Fig.  3a  clearly  suggests  that  equal  partitioning  is  not  occurring. 

Figure  4  is  the  same  as  Fig.  3,  except  the  one-way  approximation  to  Eq.  1  was  used  to  compute  the 
modal  amplitudes.  The  one-way  approximation  is  more  consistent  with  the  approximations  used  to 
derive  Eq.  3,  namely,  the  neglect  of  the  backscattered  acoustic  field.  Although  one  observes  qualitative 
agreement  of  Figs.  4a  and  4b,  equal  partitioning  of  modes  in  Fig.  4a  occurs  to  a  lesser  extent  than  in 
Fig.  4b.  One  might  make  the  tentative  conclusion  that  the  RHS  of  Eq.  3,  based  on  a  Bom-like 
approximation,  does  not  adequately  capture  the  physics  of  scattering  for  seabed  roughness  having  a 
small  correlation  length  as  compared  to  modal  interference  lengths  in  the  trapped  spectrum.  The 
relatively  smaller  correlation  lengths  can  inhibit  energy  from  being  transferred  from  the  trapped  to  the 
continuum,  which  in  turn  inhibits  equal  partitioning  of  modes. 

The  above  are  only  tentative  conclusions  and  require  additional  study  of  how  to  properly  apply  the 
master  equation  for  waveguides  that  exhibit  strong  variability. 
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FIGURE  1.  Modal  amplitudes  versus  range. 
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FIGURE  2.  A(KO  for  h=  0.5  m. 
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FIGURE  3.  Modal  power  flow  predicted  by  two-way  solution  in  (a)  LHS  of  Eq.  3  and 

(b)  RHS  of  Eq.  3. 


FIGURE  4.  Modal  power  flow  predicted  by  one-way  solution  in  (a)  LHS  of  Eq.  3  and 

(b)  RHS  ofEq.  3. 
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IMPACT/APPLICATIONS 


One  potential  impact  of  this  research  is  that  these  studies  may  assist  in  understanding  how  to  optimally 
combine  advanced  propagation  models  (non-separable  and  3-D)  and  information  inference  methods  as 
one  proceeds  to  study  ocean  waveguides  with  increasing  complexity  and  inhomogeneity. 

TRANSITIONS 

The  combined  study  of  statistical  inference  and  the  effects  of  seabed  scattering  is  expected  to  relate 
propagation  statistics  to  physical  mechanisms  on  continental  shelf  and  slope  environments.  The 
knowledge  of  such  statistics  and  their  relationship  to  the  physics  of  the  propagation  may  be  useful  for 
sonar  applications  in  continental  shelf  and  slope  environments. 

RELATED  PROJECTS 

Related  research  projects  include  applying  maximum  entropy  to  the  statistical  inference  of  values  for 
seabed  geoacoustic  and  scattering  parameters  from  measured  data. 
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